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ABSTRACT 

In this paper we show a plausible mechanism that could lead to the formation of the Dark 
Lanes in Lunar Swirls, and the electromagnetic shielding of the lunar surface that results in the 
preservation of the white colour of the lunar regolith. We present the results of a fully self-consistent 
2 and 3 dimensional particle-in-cell simulations of mini-magnetospheres that form above the lunar 
surface and show that they are consistent with the formation of ‘lunar swirls’ such as the archetypal 
formation Reiner Gamma. The simulations show how the microphysics of the deflection/shielding 
of plasma operates from a kinetic-scale cavity, and show that this interaction leads to a footprint 
with sharp features that could be the mechanism behind the generation of ‘dark lanes’. The physics 
of mini-magnetospheres is described and shown to be controlled by space-charge fields arising due 
to the magnetized electrons and unmagnetized ions. A comparison between model and observation 
is shown for a number of key plasma parameters. 


Popular Summary 


On the near side of the Moon, located on one of the dark, most ancient plateaus of the Moon, 
there is an unusual white, swirl marking that appears to bear no relation to any typical impact 
craters or lunar ejecta. Visible from Earth, the Reiner Gamma Formation, is the quintessential 
example of a lunar swirl. 

The lunar surface darkens with age due to space weathering from the bombardment of solar wind 
protons. The presence of small, lO-lOOOkm, regions of unweathered lunar surface where no topo¬ 
graphical cause can be determined, have been an intriguing astronomical phenomena for centuries. 

During the early Apollo missions it was discovered that the location and size of lunar swirls, 
coincide with small (~ 100s km) surface patches of crustal magnetic field. The Moon is without 
either an atmosphere or an overarching magnetic field, unlike the Earth’s, that would erode and 
disperse differences in rates of space weathering maturation. The lunar swirl markings therefore 
represent a time integrated record of the ability of low intensity magnetic fields (~ l/100 th that 
at the Earth’s surface) to shield parts of the Moon’s surface from proton bombardment. The fine 
detail visible within the swirls points to a very precise and consistent control of the proton flux. 
Understanding how this works could lead to the ability to artificially recreate and enhance the 
protective effects for people and instrumentation susceptible to damage from the more intense parts 
of the cosmic radiation. 

In this paper 2D and 3D particle-in-cell plasma simulations are used to show how finite Larmor 
orbit plasma effects reproduce the small scale deflection, reflection and retardation of incoming 
interplanetary plasma with the detail level of the lunar swirls and the plasma observations recorded 
by surveying spacecraft. 


Introduction 


The Moon does not, and may never have had, an ac¬ 
tive core dynamo with which to generate its own global 
magnetic field [I]. During survey missions for Apollo, 
however, it was discovered that there were several small, 


static regions of magnetic field on the surface of the 
Moon |2J [3j. The dimensions of these regions of mag¬ 
netic field are of the order of 100s of kilometres. This 
suggests that they originate from magnetised material 
on, or very close to, the lunar surface. These small 
patches of crustal magnetic field have been identified with 
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the formation of “miniature magnetosphere” cavities 
[6] and show many observational characteristics of mini¬ 
magnetosphere boundaries [M3. This is despite their 
relatively small dimensions, many orders of magnitude 
smaller in dimensions than their planetary cousins, and 
much lower magnetic field intensity (believed to be only 
of the order of 100-300nT at the surface). Particular in¬ 
terest in these features is provided by the apparent link 
between the mini-magnetospheres and the creation of ‘lu¬ 
nar swirl’ patterns on the Moon mm- 

Lunar swirls are optically distinct, white, curvilinear 
surface features that are found in several locations across 
the lunar surface HZ!- The features are distinctive due to 
their fluid or wispy structure that is unlike either craters 
or impact ejecta. Their form has been determined to 
be unrelated to geographical topography and appears to 
overprint on both mountainous and plateau terrain [UK 
[20]. What is found associated with all lunar swirls is 
that they are always co-incident with similar sized areas 
of magnetic field [20]. However not all the anomalous 
crustal magnetic fields on the Moon have identifiable lu¬ 
nar swirl discolourations [2T1 . 

Lunar swirls and magnetic anomalies 

The presence of small areas of magnetic field on an oth¬ 
erwise unmagnetised planetary body is not unique [23- 
124] . The distribution of magnetic field anomalies on the 
Moon, vary from irregular conglomerations and clusters 
to small and isolated - such as the Reiner Gamma for¬ 
mation [H[25Mj- The largest distributions of crustal 
anomalies are located on the southern part of the far- 
side of the Moon, antipodal to the Crisium, Serenitatis, 
Imbrium, and Orientale basins [HESM1 • 

Several theories exist to explain the creation of of 
swirls [23[ [28[ [29]. However, recent work •'{(] on data 
from the Lunar Reconnaissance Orbiter Diviner Lunar 
Radiometer support the hypothesis m 125 that the 
mechanism of the variations in albedo is related to dif¬ 
ferential solar wind particle bombardment of the lunar 
regolith. 

The continuous bombardment by solar wind ions and 
micrometeorites alters the surfaces of all airless bodies in 
the solar system m- The effect of this space weath¬ 
ering include an overall decrease in visible to nearIR re¬ 
flectance, attenuation of mafic absorption features and 
introduction of a strong positive slope (spectral redden¬ 
ing) [32]. These changes are attributed to the production 
amorphous coatings on grain surfaces and tiny blebs of 
metallic iron, known as nanophase iron (npFeO) [331 134| . 
The corollary of this process is that a reduced particle 
flux should lead to lighter coloured (i.e.unweathered) re¬ 
golith, and an enhanced flux should lead to even darker 
regolith than the backgound [20] [2TJ [30] El] ■ 

The fine details visible within the curvilinear shapes of 


the swirls is often accentuated by dark lanes that wind 
though the bright swirls. The transition between the 
light and dark is often very sharp and the width of the 
dark lanes can be narrow, less than a kilometre. These 
dark lanes suggest a consistent enhancement in proton 
bombardment relative to on-swirl reduction and near-by 
off-swirl surface representing ‘normal’ lunar regolith m 

El EM- 

For this level of detail to have formed and been re¬ 
tained by the swirls, points to a very precise and consis¬ 
tent control of the proton flux. The finesse of the tran¬ 
sitions further suggests the process occurs close to the 
surface. A remote plasma structure would be more likely 
to be shifted and dispersed by the fluctuations in the so¬ 
lar wind and lunar cycle plus effects of transitions in and 
out of the regions of the Earth’s magnetotail. 

The question now becomes determining the mechanism 
for the consistent, localised, low altitude redirection of 
the solar wind protons that could produce the distinct 
patterns of the swirls. This depends on the fundamental 
plasma physics that determine the electromagnetic struc¬ 
tures arising as a result of the solar wind impacting the 
small scale magnetic anomalies. 

The characteristics of this interaction above the sur¬ 
face have been measured by several in-situ spacecraft [41 
El El EM- This allows direct comparisons to be made 
between models and measurements. It is clear from 
the data that the solar wind ions are reflected and de¬ 
flected by the magnetic fields as with planetary mag¬ 
netospheres [4]. Conclusive evidence that the crustal 
magnetic anomalies are also capable of producing colli¬ 
sionless shocks has also recently been demonstrated by 
the analysis of Halekas and co-workers Ell¬ 
in determining the overall form of a planetary mag¬ 
netosphere, the ions and electrons of the surrounding 
plasma environment can be approximated as being locked 
together with each other and the magnetic field lines. 
However, the small size of the magnetic ‘bubble’ cre¬ 
ated by the lunar magnetic anomalies is larger than the 
scale of the electron dynamics (< 1 km— ~20 km), but 
very much smaller than the ion dynamics scale (~100s- 
1,000s km). This makes lunar magnetic anomalies ideal 
natural laboratories for studying fundamental collision¬ 
less plasma physics phenomena and to validate the small 
scale kinetic processes in particular, finite ion Larmor 
radius effects of magnetised, collisionless shocks. The 
equivalent sizes in magnetic confinement or inertia fusion 
plasmas would be ~ 1mm and ~ 1 /im. 

In a previous paper ;30], in-situ satellite data, theory 
and laboratory validation showed that it is an electric 
field associated with the small scale collisionless shock 
that is responsible for reflecting, slowing and deflecting 
the incoming solar wind around mini-magnetospheres. It 
was shown that the electric field of polarization, caused 
by the gradient in the magnetic field, between charge car¬ 
riers of the solar wind flow, is of prime importance. This 
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polarization field leads to reflection and scattering of the 
protons and electrons [40]. The counter-streaming of 
the ions ahead of the barrier is responsible for generating 
lower-hybrid waves via the modified two-stream instabil¬ 
ity ED, a kinetic plasma process. 

In this paper our 2- and 3D simulations support the 
observational evidence and analytical arguments this re¬ 
sults in a self consistent explanation for the creation of 
‘lunar swirls”. It will be shown how the 3D particle-in¬ 
cell plasma simulation reproduces the key observational 
characteristics from both above and on the surface, using 
the simplest of magnetic tomography models - a single 
magnetic dipole. Thus illustrating the complexity ob¬ 
served comes from the plasma-magnetic field interaction 
at these scale sizes. In the 2D simulations shown, some 
observables - such as barrier width - can be seen to be 
consistent irrespective of changes to magnetic field ori¬ 
entation, size and solar wind conditions. Conversely the 
variability of other characteristics - such as stand-off dis¬ 
tance - can alter significantly with environmental plasma 
conditions as expected. The extended swirl pattern can 
be envisaged as being composed from a scattered distri¬ 
bution of the elemental components (a subset of which 
are shown here as examples) with varying extent, orien¬ 
tation, field intensity and overlap. 

Other features, such as the most distinctive “eye” of 
the Reiner Gamma formation [42], will be shown to 
be consistent with a single dipole (in 3D), laid on it’s 
side with it’s magnetic axis parallel to the lunar surface, 
interacting with a downward flowing solar wind plasma 
in accordance with the theoretical principles outlined in 
[39] and computationally illustrated here. A recreation 
is shown of exaggerated typical spacecraft instrumen¬ 
tation signatures for a transit through a simple mini¬ 
magnetosphere. 

METHOD 

The plasma simulation was carried out using a particle- 
in-cell (PIC) code, called OSIRIS [15] . In OSIRIS the 
full set of Maxwells equations are solved on a grid using 
currents and charge densities calculated by weighting dis¬ 
crete particles onto the grid. Each particle is pushed to a 
new position and momentum via self-consistently calcu¬ 
lated fields. The code makes few physics approximations 
and is ideally suited for studying complex systems with 
many degrees of freedom such as this one. The cost of this 
fidelity is the computational expense of the larger num¬ 
ber of more complex algorithms interacting with each 
other than in a simulation of a more simple model sys¬ 
tem. The reason this is necessary is that the scale size of 
the mini-magnetosphere structure is very much smaller 
than the hydro-code approximations allow. No filtering is 
performed in these simulations, this allows us to resolve 
all waves in space and time including whistler waves. The 


code is a time and space domain code, not a spectral code, 
so the equations are integrated via FFTs. 

Mostly shown here are the 3D simulations with cut- 
throughs (see graphic in Figure [l]) to highlight some fea¬ 
tures. The full 3D simulation is important because plas¬ 
mas inherently have orthogonal acting forces and flows 
such as E = v x B, and some aspects are missed in 2D 
simulations (examples of which are shown in Figurefj7]to 
|)]). When the dimensions of the magnetic anomalies are 
compared to the characteristic scale length of the mag¬ 
netic inhomogeneity/magnetic anomaly, Xb , then the 
value obtained for Xb = [(1 /B sw ).(d,B/dh)]~ l is ~ lkm, 
this is assuming a surface magnetic field of 200nT and 
the boundary to be at an altitude, h, of 20km. 

What this means is the incoming ions have gyro radii 
of the order of the structure size. Therefore their motions 
are not simply following the magnetic field lines, they are 
unmagnetized. However the electrons gyro radius is suf¬ 
ficiently small, compared to the overall size of the struc¬ 
ture, that they are able to follow the changes in magnetic 
field created by the magnetic anomaly. The difference in 
behaviour between the electrons and ions sets up a space 
charge that controls the ions behaviour [44] (45i . To rep¬ 
resent the dynamics it is necessary to perform either lab¬ 
oratory simulations |39i |3Hj or full PIC code simulations, 
like OSIRIS, as is done here. 

The simulation code operates in normalised plasma 
units. The independent variable being the plasma den¬ 
sity n sw . The units for distance are normalised multi¬ 
ples of the electron skin depth, c/oj pe (where c is the 
speed of light and u> pe is the electron plasma frequency 
u) pe = \Jn e e 2 / m e e 0 and e the electric charge, m e the elec¬ 
tron mass and e 0 the permittivity of free space). From 
which it can be seen that the physical dimensions of the 
box scale inversely with the square root of the free pa¬ 
rameter of plasma density, n sw . The equivalent dimen¬ 
sions in kilometres recreated by the simulation can be 
made to contract with higher densities and expand with 
lower as other de-normalised parameters, such as mag¬ 
netic field intensity, would also increase and decrease to 
remain in proportion with the defined important plasma 
parameters (specifically the magnetosonic mach number 
and plasma /3) that determine the correct realm of plasma 
environment on the Moon. 

Simulation parameters 

The 3D simulation geometry and initial conditions are 
shown in Figure [I] The lunar surface is represented by 
the lower x — z plane. A single source magnetic dipole is 
placed just below and parallel to the x — z plane, with 
magnetic axis aligned along x = 0, with the north pole 
orientated in the +z direction. The result is a hemi¬ 
spherical magnetic field emerging from the surface. A 
magnetised ‘solar wind’ plasma with density n sw and 
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magnetic field B sw (antiparallel to anomaly field) is in¬ 
troduced from the top plane with a flow velocity —v sw 
vertically down onto the lunar surface. 

Observationally the mean value of the solar wind mag- 
netosonic Mach number is M m = 8 [47], with a /3 = 0.2, 
n sw = 10cm -3 , B sw = lOnT and Ti = 5eV. The electron 
and ion skin depths are = 1.7km and 72km respectively. 
The ion Larmor orbits are 97 km for the thermal distri¬ 
bution and 627km for the flow velocity. 

This makes c/uj pe = 1.7 km and the simulation box 
shown represents ~ 340 x 340 km by ~ 170 km altitude. 
The grid cell size is 5 cells electron skin depth, i.e., grid 
resolution of ~300 m. In this simulation the magnetic 
moment of buried dipole = 9.10 15 A m 2 at a depth of 
~50 km (for density of n sw = 10 cm -3 ). 

The tradeoff for preforming full-pic in 3D required op¬ 
erating within the simulation the plasma flow speed is 
increased by a factor F = v sw /v OS i r i S . In the 3D simula¬ 
tions F = 20 assuming a v sw = 600 km s -1 . 

In order to maintain dynamic similarity with the lu¬ 
nar environment, the magnetic fields and temperatures 
were scaled proportionately so as to maintain the same 
control variables of plasma /3 (thermal pressure to mag¬ 
netic pressure) and M m (magnetosonic mach numbers 
and speeds). The fast magnetosonic mach number M m = 
v sw /\/ v \ + v s where v s is the sound speed and va is 
the Alfven speed. For the same plasma density n, the 
temperature within the simulation was T OS i r i a = T SW F 2 
and magnetic field intensities was B 0 siris = B SW F. The 
plasma frequency cj pe = \Jne 2 /me Q is unaffected. The 
real proton-to-electron mass ratio of 1836 was used. 

For this simulation the peak magnetic field at the sur¬ 
face, B s f, scales as ~ 2000nT for a realistic quiet sun 
solar wind of v sw = 600 kms -1 and n sw = 10cm -3 . 

The standoff position r s , in the 3D simulation, of the 
cavity was engineered to be at the ion skin depth c/uj p i 
(for a plasma density of 10c?n 3 , r s = 100km, B sur f ace ss 
2/rT). The motivation for this was that observation of the 
magnetic compression would occur at the same altitude 
and completely within the simulation box. 

The consequences of the scaled up velocity are that 
the simulation surface magnetic field had to be similarly 
scaled up to match the flow pressure. Also not all the 
microphysics scales, in particular excited waves, collision¬ 
less dissipation mechanisms may be reproduced. But the 
microphysics of the collective space-charge electric held, 
which is responsible for the plasma shielding/deflection 
is well captured. This then results in a sharp foot print 
with features on ~ c/w pe scales as will be shown later. 

The analytical expression for r s related to the pres¬ 
sure balance between incoming solar wind n sw mv 2 w = 
B 2 / 2/u 0 , is well established for planetary e.g. [53J, and 
mini-magnetospheres [48]. Such that: 



Here B a is the subselenean dipole moment. 

Using realistic magnetic held strengths will lead to 
smaller cavity structures and lower r s . 

SIMULATIONS RESULTS 

The simulation code is used here to recreate simpli¬ 
fied case study combinations of plasma parameter condi¬ 
tions and magnetic held dipole orientations and intensi¬ 
ties. These can be compared to the analytical expressions 
[55] and observational data e.g. m- Included here are the 
importance of simulating 3 dimensions in a plasma due 
to the inherent cross-product relationship of forces. 

Figure [2] shows the results of the simulation of the so¬ 
lar wind plasma impacting a localised crustal magnetic 
held structure. In Figure [2] only the ion density above 
a threshold is visualised in order to make the box trans¬ 
parent. The magnetic held structure is shown by the 
blue held lines (again in-part omitted for sake of clar¬ 
ity). The red represents the space-charge electric held at 
the boundary, which is set up by the different penetra¬ 
tion depth between ions and electrons at the edge of the 
magnetosphere. The lateral-projections of the electric 
held structure (corresponding also to the relative proton 
density) reveal interesting dynamic features and orthog¬ 
onal asymmetries. The projection on the y — z plane 
shows a rippled surface structure, due to the diamagnetic 
electron-ion drift instability which occurs perpendicular 
to the magnetic held lines [55] ■ In contrast, the projec¬ 
tion on the x — z plane shows a smooth surface structure 
since the relative electron-ion drift is absent. This il¬ 
lustrates the anisotropic preferences of particular plasma 
instabilities. However the narrow width of the barrier 
remains a consistent feature, although not necessarily a 
single, smooth boundary due to waves, turbulence and 
instabilities - the magnitude of which alter with specihc 
conditions. 

Representative trajectories of a few solar wind protons 
are shown by the small spheres and yellow track lines 
and are seen to be widely scattered not by a gradual 
redirection but ballistically scattered from a very narrow 
region close to the surface. The electric held responsible 
for scattering can be seen to be omni-directionally point¬ 
ing outwards regardless of the magnetic held orientation. 
This is because it relates to the gradient in the magnetic 
held intensity V|U| 2 and not B as is predicted theoreti¬ 
cally [39]. The projection onto the x — z or surface plane, 
shows the electric held intensity at the lunar surface. The 
proton density is controlled, on these scales, by the elec¬ 
tric held rather than the magnetic held because the ions 
are unmagnetised. The simulation here, therefore, shows 
the distribution of sharp regions of enhanced proton-flux 
and regions of depleted proton hux. 

In order to reveal the details of the interior structures 
of the 3D simulation Figure [3] shows the 2 dimensional 
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sections of each of the plasma parameters. These are 
from top to bottom the ion and electron density m, n e , 
the resultant electric E and magnetic field B and the 
original source half dipole magnetic field undisturbed by 
the solar wind plasma B vac . The plane of the 2D sec¬ 
tion bisects the midpoint in x of the dipole axis and is 
indicated in Figure [l]. 

The density pile-ups, exclusion of the majority of the 
particles from the interior, back-flow and turbulence of 
the barrier are apparent on all parameters. The barrier 
currents that form result in further responsive inductive 
currents and corresponding inductive magnetic fields. 
The thickness of the barrier is of the order of the elec¬ 
tron skin depth ~ c/w pe , as theoretically predicted TfDl . 
The small-scale plasma instabilities, waves and turbu¬ 
lence formed, provide the means by which the ion and 
electron particle distributions become non-thermal and 
exchange energy through Landau damping [451 . (This is 
illustrated in the top two panels of Figure [lO] discussed 
in Section I. 

The “stand-off” distance, r s (equation]]) that the mag¬ 
netic cavity reaches force balance with the incoming 
plasma, is approximately 50 c/oj pe by design so as to en¬ 
sure the interface was well situated within the simulation 
box. 

The relative proton density from the simulation can be 
further clarified by examining a linear plot taken at the 
equivalent lunar surface level. This is shown in Figure [4] 

If sustained long term, this pattern of excluded and 
narrow enhance proton density resulting from the action 
of the mini-magnetosphere with crustal magnetic field 
would be consistent with the distribution and finesse re¬ 
quired to form the variety of lighter and darker albedo 
alterations seen in lunar swirls. 

The relative deposition of proton flux (shown in grey) 
on the surface slice from the simulation is shown in the 
top left image of Figure [5] with the magnetic field lines 
of the dipole. This is compared to a image of the central 
region of the most distinctive example of a lunar swirl, 
Reiner Gamma Formation located at 7.4°TV, 300.9° E, 
taken by NASA’s Lunar Reconnaissance Orbiter 50]. 
The agreement between the key characteristics of dark 
lane [5L width and shielded interior can be seen to be 
totally consistent. The symmetries in the simulation re¬ 
sults prevent distinguishing categorically the orientation 
of the magnetic fields. 

We will now preform a series of simulations to study 
the parameter regime of single dipoles by varying the 
orientation and some of the plasma conditions. 

The effect of dipole orientation 

Changing the orientation of the emerging magnetic 
dipole relative to the surface plane produces the same 
diamagnetic characteristics of narrow barrier, particle re¬ 


flection, cavity formation, waves and turbulence. How¬ 
ever, as can be seen in Figure [6] the overall morphologies 
of the relative proton distribution are very different. 

Figure [6] shows 2D simulations of the same magnetic 
dipole in 3 orthogonal directions x, y and z, relative to 
the surface plane. The solar wind, mass ratio and plasma 
conditions are as those of the 3D simulation with the 
exception that the solar wind plasma is not magnetised. 
The consequence of B sw parallel or antiparallel to the 
surface magnetic field is shown in the next section. 

The linear plots of relative proton density (on the right 
of the Figure) show different proton deposition profiles 
for each orientation. 

The affects of changes in the environmental plasma 

Although the crustal magnetic anomalies are a fixed 
magnetic field source, the plasma environment is not. 
There are periodicities due to orbits and diverse solar 
wind and/or magnetospheric conditions. Through simu¬ 
lation this can be explored, albeit briefly, in Figures [7] [H] 
and [H 

A spacecraft passing over the vicinity of these anoma¬ 
lies, at a fixed height (as indicated in the Figures by a 
satellite graphic and dashed ‘flight-path’ line) but under 
the different plasma conditions will sometimes transect 
the different regions of a mini-magnetosphere and so ob¬ 
serve different characteristics. (An illustration of these 
plasma data instrumentation signatures in such a flyover 
is shown in Section [|. 

Figures [7] and [8] illustrate the consequence of combi¬ 
nations of changes the simulation conditions. In Fig¬ 
ure [7] the incoming solar wind Magnetosonic mach num¬ 
bers M m and parallel (or anti-parallel) solar wind mag¬ 
netic field B sw orientation are varied. In Figure [8] the 
dipole size/length, Ldipole, is altered. The extent of the 
dipole can be altered either through being a larger single 
crustal magnetic field or by multiple smaller conglomer¬ 
ations (such as on the far-side of the Moon) that will 
appear as a single dipole when observed in sufficient al¬ 
titude to be in the far-held. 

In all these figures the white-to-red colour distribution 
represents relative ion density and all the conditions are 
the same except for those stated as otherwise. 

The majority of the simulation results shown have been 
with the solar wind plasma impacting normal to the plane 
of the lunar surface andemerging magnetic held. This is 
purely for simplicity and consistency so as to example 
the plasma interaction. In Figure[9]the impacting plasma 
how is parallel to the lunar surface and so can be seen to 
create a mini-magnetotail. From Figure [9] it can be seen 
that many of the features of the normal instance impacts 
are also seen in the parallel impact case. However the 
diversity of regimes and features that can be discussed 
is greatly increased and more suitable for a dedicated 
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study. This is also true of combinations of magnetic field 
anomalies and conditions. 

In summary, the 2D simulations show how only occa¬ 
sionally might certain features, like the diamagnetic cav¬ 
ity, be clearly detected in in-situ instrumentation, such 
as particle detectors and magnetometers. This is espe¬ 
cially for a rapid spacecraft transit through very small 
isolated magnetic held anomalies. However, other instru¬ 
mentations such as imagers [5], that look down onto the 
features could still detect the characteristic depletion in 
proton reflection. 


Simulating Spacecraft flyover signatures 


Figure [To] shows how the simulation results translate 
as observations to a spacecraft flying above the surface 
anomaly. 

These features should be common to any mini¬ 
magnetosphere to some degree or another irrespective of 
the presence of lunar swirls. The comparison is qualita¬ 
tive not quantitive to highlight the nature of interplay 
in the parameters so as to provide identification of phe¬ 
nomena between theory and observation via simulation 
visualisation. 


The plasma signatures reproduced in the simulation 
are: (a) An increase and density pile-up before the en¬ 
try into the cavity (reported observationally by [51 37, 
[551 152)). This is most evident in Figure 10 panels la¬ 
beled ion (m) and electron (n e ) densities, (b) The mag¬ 
netic field components (panels labelled B z and B x in Fig¬ 
ure [TO]) are seen to rotate as if “draped” around a small 


magnetic obstacle as first reported by Lunar Prospec¬ 
tor [52] . (c) An increase in electrostatic (E) held and 
solitary waves at lower hybrid frequencies appear (obser¬ 
vational data [6® 12 [53]). (d) Non-adiabatic energy be¬ 
tween the ions and electrons (37], [38]. In the upper most 
panel, the electron momentum p z e in the vertical/how 
(z) direction, increases and the energy distribution (tem¬ 
perature) changes as seen by the spreading in p Zi6 . (e) 
Back streaming protons (37} [38] accelerated by similar 
factors close to the shock surface [6]. The protons mo¬ 
mentum p z i (next panel) shows how the ions are reflected 
back from the magnetic structure predominantly at in¬ 
coming velocity, (f) Electrostatic oscillation at the lower 
hybrid frequency (seen observationally |55j). The back- 
flow of ions establishes the conditions for the modified 
two stream instability that oscillates at the lower hybrid 
frequency m- 


The simulation operated with a v sw x 100 the typical 
realistic quite time velocity of 600 kms -1 and a reduced 
mass of 100. This emphasised in the simulation, the non- 
adiabatic energy between the ions and electrons (top two 
panels). 


CONCLUSIONS 


We have reproduced all the major the characteristics 
of lunar swirls using the simplest of magnetic topologies 
- a single dipole. 

The simulations confirm many of the satellite find¬ 
ings and the theoretical predictions [5] that a collision- 
less shock forms at the altitude expected from the the¬ 
ory of collisionless shocks. The thickness of this shock 
is approximately equal to the electron skin depth c/w pe 
(where u> pe is the electron plasma frequency) again in 
agreement with theory. The characteristic observations 
of electron and ion density enhancements and depletions 
accompanied by magnetic field intensity pile-up at the 
shock coincide with the formation of a narrow interface 
region where a dynamically stable electric held exists. 
The simulations confirm that it is this electric held that 
controls the behaviour of the solar wind ions impacting 
the magnetic structure. 

Because the primary driving term is a gradient in en¬ 
ergy density, it is not exclusive to one magnetic held ori¬ 
entation. There are asymmetries in the diamagnetic cav¬ 
ity due to a difference in the preferred plasma instabilities 
and different growth and instability interchange rates. 

The 3D simulation has the solar wind howing from di¬ 
rectly above the anomaly, with normal incidence to the 
surface. The interaction naturally leads to an asymmet¬ 
ric structure in the helds and density because the plasma 
electrons are preferentially deflected to one side of the 
magnetic held rather than the other. Though the mag¬ 
netic structure is symmetric, the plasma response is in¬ 
trinsically asymmetric. This does not include a magneto¬ 
tail from an oblique impact, however it is the most basic 
prototype example. 

In the case of the Moon, different orientations of the 
magnetic held structure to the (a) surface of the Moon 
and (b) to the incident solar wind - will result in more 
complex and filamentary structures, particularly near 
magnetic poles. The size of the hlamentary structures 
could be less than 1 km or even 100 m. This allows a 
complex surface magnetic held to provide variable levels 
of protection from the ageing process of the solar wind. 
What is seen on the surface is the legacy of structures 
above and the periodicities and variations encountered 
over time. 

A prediction from this work is that the dark lanes 
should be always approximately the same width (~ 
c/wp e ) everywhere on the Moon. Whereas the extent 
of the white regions will depend upon magnetic held 
strengths and near held topography. Future work is to 
determine the near-held tomography, ideally with in-situ 
measurements on the ground. We will also look at more 
complex magnetic configurations and combinations plus 
mini-magnetotail formation in a future papers. 

Although not all the microphysical phenomena (e.g. 
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dissipation, excited waves, etc) scales correctly in the 
simulation nor rely on these simulations to extract re¬ 
alistic particle distribution data, and exact growth rates 
of the excited waves upstream of the reflection of the 
plasma, etc. Although these are first principle simu¬ 
lations, they are still unable to capture all the physics 
quantifiably exact, because it is still currently impossible 
to simulate the realistic physical parameters due to the 
massive computational requirements. 

These simulations do illustrate clearly, however, the 
microphysics underlying the shielding/deflection of the 
plasma around a kinetic-scale cavity. 
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FIG. 1: The initial conditions of the simulation. The magnetic field intensity is shown projected onto the back walls 
and ground plane in banded grey where decreasing band interval corresponds to increasing magnetic field intensity. 
Selected magnetic field lines are shown in graduated blue. The magnetic dipole moment m m is 25 normalized units 
long and resides 25 units below the center of the box aligned with the positive x axis. The plane marked ‘2D slice’ 
shows the relative orientation of the sections shown in Figures [3] and [4j 



FIG. 2: (Colour) A 3D magnetized plasma collision with a surface magnetic dipole. The solar wind plasma (flowing 
vertically downwards) impacts a localised crustal magnetic field structure (blue lines). The green spheres and tracks 
show a subset of protons population trajectory being scattered from the narrow polarisation electric field (red). 
Proton density enhancements follow the electric field. Only part of the magnetic field lines are shown for clarity and 
the background densities are not visualised. 
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FIG. 3: (Colour) Planar slices through the 3D simulation shown in Figure [ 2 ] showing the behaviour of the plasma 
densities and electromagnetic forces. The location of the plane is indicated in Figure |T) The normalised y axis 
represents altitude above the Moons surface. The normalised x axis represents distance along the surface of the 
Moon upon which a magnetic dipole field is located. 
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FIG. 4: The simulated relative proton deposition onto the surface of the Moon. From the same simulation shown 
Figures 1|2|3 



FIG. 5: Left:A 2D slice of the relative proton density from the 3D simulation with the initial magnetic field lines 
from a single subsurface dipole. The greyscale distribution shows darker for higher density of protons, whiter for 
less. Right: A image of the central region of the Reiner Gamma Formation lunar swirl taken by NASA’s Lunar 
Reconnaissance Orbiter [50].The form and relative width of the ’dark lanes ’ m suggest the aspect ratio of 
dark-lane width to cavity width is similar in both cases. 
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FIG. 6: Different orientations of magnetic dipole axis aligned in the x,y,z direction with their corresponding relative 
proton density “footprint on the lunar surface. The double high peak visible in the vertical dipole orientation 
(middle panel) would in 3d, be the equivalent of the particle deposition pattern that would form an auroral oval. 



FIG. 7: The affect of changes in parameters 1. Perpendicular shocks by cold plasma with fixed dipole size/strength 
(Ldipoie = 2), varying magnetosonic Mach number, M m = 3,5, for parallel and anti-parallel Bdi po i e • Reduced 2d 
parameters: The incoming flow velocity is x50 a typical solar wind speed, proton to electron mass ratio m p /m e is 
~ 1/20 realistic value, chosen for computational speed. Although these values are not representative of actual values 
in space, they do allow a qualitative comparison of the variation of mini-magnetosphere characteristics. 






























































13 


B sw antiparallel 



-50 0 50 100 

X [c/COp] © e,p Alves, R.A. Bamford 


FIG. 8: The affect of changes in parameters 2. Perpendicular shocks by cold plasma with fixed magnetosonic Mach 
number M m = 2.0 with Bi n parallel to Bdi po i e > varying dipole size/strength Ldi po ie = 0.5,2.5, 5.. Reduced 2d 
parameters as . 
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FIG. 9: The formation of a mini-magnetotail. The solar wind flow direction is parallel to the lunar surface and the 
dipole vector is pointing outwards from the page. In the upper panel the magnetic field intensity is shown as banded 
grey. 
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FIG. 10: (Colour) The simulated spacecraft diagnostic signatures of a transit through a mini-magnetosphere. The 
2D OSIRIS simulation is plotted so as to recreate the type of signatures that would be observed by spacecraft 
plasma instrumentation during a flyover of a crustal magnetic anomaly (lower most panel). A low altitude lunar 
spacecraft would record in a flyover transit over the surface anomaly at constant altitude of h = 5.5 x c/u> pe in 
normalised units the equivalent of ^12km (for a 5cm -3 density plasma). The simulation is the simplest geometry 
with slow wind flow normal to the lunar surface, there is no drawn out magnetotail in this example. The simulated 
data window would be the equivalent of ~2 to 4 minutes in duration. 


































